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ABSTRACT 

We  prove  an  existence  theorem  locally  in  time  for  quasilinear  hyperbolic 
equations,  in  which  the  coefficients  are  allowed  to  depend  on  the  history  of 
the  dependent  variable.  Singular  perturbations,  which  change  the  type  of  the 
equation  to  parabolic,  are  included,  and  continuous  dependence  of  the  solu¬ 
tions  on  the  perturbation  parameter  is  shown.  It  is  demonstrated  that,  for  a 
substantial  number  of  constitutive  models  suggested  in  the  literature,  the 
stretching  of  filaments  of  polymeric  liquids  is  described  by  equations  of  the 
kind  under  study  here. 
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SIGNIFICANCE  AND  EXPLANATION 


In  an  earlier  report  [23],  I  proved  an  existence  theorem  locally  in  time 
for  a  class  of  parabolic  differential-delay  equations  modelling  the  stretching 
of  filaments  of  polymeric  liquids.  For  many  materials,  however,  the  term 
determining  the  parabolic  nature  of  the  equation  seems  to  be  very  small.  The 
question  was  therefore  raised  whether  an  existence  theorem  could  still  be 
proved  in  the  limit  when  this  term  tends  to  zero.  This  paper  gives  a  partial 
affirmative  answer  to  that  question.  For  a  broad  class  of  constitutive  assump¬ 
tions  suggested  in  the  rheological  literature,  the  problem  can  be  described  by 
quasilinear  hyperbolic  equations.  We  prove  an  existence  theorem  for  these, 
and  we  also  prove  continuous  dependence  on  small  parabolic  perturbations.  The 
results  are  based  on  theorems  due  to  Kato.  In  one  respect  the  problem  studied 
here  is  different  from  the  one  in  [23] s  Instead  of  a  filament  pulled  at  its 
ends,  we  study  an  infinite  filament  under  the  influence  of  a  longitudinal  body 
force.  It  is  hoped  that  similar  results  can  eventually  be  obtained  for  the 
former  -  physically  more  relevant  -  problem  and  also  for  more  general  flow 
geometries.  _ _ 
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The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


SINGULARLY  PERTURBED  HYPERBOLIC  EVOLUTION  PROBLEMS  WITH  INFINITE  DELAY 
AND  AN  APPLICATION  TO  POLYMER  RHEOLOGY 
Michael  Renardy 

1 .  introduction . 

In  a  recent  paper  [23],  I  proved  an  existence  theorem  (locally  in  time) 
for  solutions  to  a  class  of  quasilinear  parabolic  differential-delay  equations 
that  can  be  used  to  model  the  stretching  of  filaments  of  polymeric  liquids. 

Such  equations  arise,  if  the  constitutive  law  is  such  that,  besides  an  "elastic" 
part  which  is  a  functional  of  the  strain  history,  the  stress  has  also  a 
Newtonian  part.  For  many  materials,  e.g.  molten  polyethylene,  however,  this 
latter  contribution  is  small.  This  warrants  a  theory  that  can  treat  the 
Newtonian  contribution  as  a  perturbation  rather  than  as  the  "leading”  term  in 
the  equation. 

In  the  present  paper,  I  shall  give  a  partial  solution  to  this  problem. 
Mathematically,  we  are  concerned  with  differential  equations  of  hyperbolic  type 
with  a  small  perturbation  changing  the  type  to  parabolic.  A  mathematical  theory 
applicable  to  such  problems  was  developed  by  Kato  [12],  [14-16].  (My  results 
in  [23]  were  based  on  the  theory  of  Sobolevskii  [25]).  Since  Kato's  theory  is 
more  easily  applicable  to  pure  Cauchy  problems  than  for  mixed  initial-boundary 
value  problems  (some  results  concerning  the  latter  are  in  [16]),  we  confine 
our  attention  to  the  former  class  of  problems  here.  Physically,  this  means 
that  rather  than  a  filament  pulled  at  its  ends,  we  will  study  the  deformation 
of  infinite  filaments  subjected  to  longitudinal  body  forces.  It  is  hoped  that 
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further  research  will  lead  to  similar  results  for  the  boundary  value  problem 
and  also  for  more  general  (in  particular  more  than  one-dimensional)  flow  geome 
tries.  It  is  clear  that  the  results  we  obtain  apply  to  other  one-dimensional 
problems  in  continuum  mechanics,  e.g.,  those  discussed  in  [6],  [9]. 

In  Chapter  2,  I  quote  those  results  of  Kato's  theory  that  are  needed  in 

this  paper.  One  of  Kato's  results  will  be  mildly  generalized.  In  Chapter  3, 

these  results  are  applied  to  a  class  of  singularly  perturbed  quasilinear  hyper 

bolic  differential-delay  equations,  which  have  the  following  form 

(1.1)  Pu  =  n  •  (f  •  u  ^  1  +  h  •  u  +  k  +  <J>,  x6R. 

tt***t  xt***t'x  xxt^'t 


n  n-1  n-2 

Here  n  is  a  small  non-negative  constant,  $  is  a  given  function  of  x 

and  t,  and  f,h,k  are  functionals  of  the  histories  of  derivatives  of  u, 

which  are  of  lower  order  than  those  displayed.  It  is  assumed  that  f  and 

h  take  posi-tive  values.  Under  appropriate  assumptions,  we  prove  that  the 

initial  history  problem  associated  to  (1.1)  has  a  unique  solution  locally  in 

time,  and,  more-over,  that  this  solution  depends  continuously  on  n, 

including  the  limit  n  ♦  0.  In  representing  differential-delay  equations  as 

abstract  evolution  problems,  we  follow  the  method  outlined  in  [23]  rather  than 

the  classical  approach  [11]. 

Chapter  4  deals  with  the  problem  of  stretching  filaments  of  polymeric 
liquids.  We  use  a  one- dimensional  approximation  to  this  problem  based  on  the 
thinness  of  the  filament  [23] .  Various  constitutive  models  suggested  in  the 
rheological  literature  [1-5],  [8],  [10],  [13],  [17-22],  [27]  are  discussed. 

It  is  shown  that,  for  all  these  models,  an  equation  of  the  form  (1.1)  is 
obtained. 

The  diversity  of  the  models  studied  here  illustrates  the  fact  that  -  at 
least  yet  -  there  is  no  particular  constitutive  law  successfully  describing 


all  the  phenomena  in  polymer  rheology.  Whether  or  not  one  constitutive  law 


can  fully  describe  a  substantial  number  of  materials,  is  not  yet  known.  As 
pointed  out  in  [26],  a  "theory  of  theories"  is  needed.  The  approach  towards 
this  goal  taken  here,  however  differs  from  that  in  [26].  Whereas  there  the 
main  emphasis  is  laid  on  results  obtainable  from  general  principles  of  con¬ 
tinuum  mechanics  or  thermodynamics,  I  found  it  worthwhile  to  look  at  the 
models  suggested  by  rheologists,  and  Z  tried  to  work  out  common  mathematical 
structures.  If  these  exist,  they  can  be  hoped  to  persist  also  in  a  more 
accurate  description,  to  which  the  models  are  more  or  less  successful 
approximations. 

One  crucial  feature  that  we  have  found  common  to  all  the  models  under 
study  here  is  that,  under  appropriate  assumptions  on  the  kernels  of  the  memory 
integrals,  it  is  possible  to  rewrite  the  equations  in  such  a  form  that  the 
highest  order  derivatives  occur  only  by  their  present  values  (in  this  context, 
cf.  [24]).  More  specifically,  for  some  n  e  H,  the  nth  time  derivative  of 
the  stress  depends  linearly  on  the  nth  and  (n+1)st  time  derivatives  of  the 
strain  with  coefficients  depending  on  the  histories  only  of  lower  order 
derivatives.  One  might  regard  this  as  a  very  general  formulation  of  the 
ancient  idea  that  polymers  combine  "elastic"  and  "viscous"  effects.  The 
oldest  models  suggested  a  linear  superposition.  What  we  have  here  may  be 
called  a  "quasi linear"  superposition. 

It  is  essential  in  our  development  that  the  integral  kernels  occuring  in 
the  constitutive  equation  are  sufficiently  smooth  everywhere,  in  particular, 
that  they  are  bounded.  This  assumption  has  also  been  made  by  other  authors 
[9],  [30],  [31].  Both  molecular  theories  and  experiments  suggest,  however, 
that  the  kernels  may  have  a  singularity  (see  the  remark  at  then  end  of  Ch. 


4).  Further  research  needs  to  be  done  on  this  point 


Acknowledgements; 

I  have  greatly  profited  from  the  lectures  of  M.  G.  Crandall,  by  which  1 
was  introduced  to  Kato's  theory.  I  also  wish  to  acknowledge  several  discus¬ 
sions  at  MRC,  which  greatly  contributed  to  motivating  this  research.  In 
particular,  I  thank  R.  B.  Bird,  C.  F.  Curtiss,  J.  D.  Perry,  M.  E.  Gurtin,  M. 
W.  Johnson,  A.  S.  Lodge,  J.  Meissner  and  J.  A.  Nohel  for  participating  in 
these  discussions. 


-4- 


2.  Abstract  Hyperbolic  Equations. 

In  this  chapter ,  I  summarize  the  results  from  Kato's  theory  that  will  be 
needed  in  the  following.  One  of  Kato's  theorems  will  be  generalized. 

We  study  an  evolution  problem  of  the  form 
(2.1)  u  »  A(t,u)u  +  f(t,u),  0  <  t  <  T,  u(0)  =  4>, 

where  u  takes  values  in  a  Banach  space  X  and  A(t,u)  is  a  linear  operator 
depending  on  t  and  u.  Our  assumptions  will  involve  further  Banach  spaces 
V  and  Z  such  that  y  c  z  c  x  with  continuous  and  dense  embeddings.  It  is 
assumed  that  V,  Z  and  X  are  reflexive  and  separable.  Let  W  denote  an 
open  set  in  Y . 

First,  we  quote  Theorem  I  from  (121  in  a  simplified  form  (with  assump¬ 
tion  N  being  obsolete)  It  is  assumed  that  the  following  estimates  hold  for 
t,t',***  6  l 0,T]  and  w,w',***  e  W  (K  denotes  a  generic  constant  indepen¬ 
dent  of  t  and  w)s 

(SI)  There  is  an  isomorphism  S(t,w):  Y  ♦  X  satisfying 

IS(t,w)»  <  K,  »S_1(t,w)Iv  <  K 

Y  rX  X , Y 

I  S(t'  ,w‘  )  -  S(t,w)ly  x  <  K(  |  t-t’  |  +  Hw-w’^) 

(A1 )  A(t,w)  generates  a  quasi-contraction  semigroup  in  X  uniformly 

„  .  A(t,w)T,  .  KT 

with  respect  to  t,w:  le  *v  *  e 

A|  A 

(A2)  S(t,w)  A(t,w)  S-1(t,w)  =  A(t,w)  +  B(t,w) 

where  B(t,w)  e  B(x),  lB(t,w)lx  x  <  K 
(A3)  A(t,w)  6  B(Y,Z)  with  «A(t,w)lv  _  <  K 

1,4. 

and  IA(t,wf)  -  A(t,w)l  <  Klwf-wl  • 

I  ,  X  X 

The  mapping  t  K  A(t,w)  6  B(y,x)  is  norm-continuous. 

(A4)  There  is  some  yQ  6  W  such  that 

A(t,w)yQ  6  Y  and  iA(t,w)y0ly  <  K 

( f 1 )  f(t,w)  6  Y,  »f(t,w)lY  <  K,  lf(t,w')  -  f(t,w)»x  <  Kl w' -wl x» 
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Moreover,  the  mapping  t  K  f(t,w)  6  X  is  continuous. 

Theorem  I  in  [12]  reads  as  follows. 

Theorem  2. 1 : 

Let  (S),  (A1)-(A4)  and  (fl )  hold.  Then  there  is  a  positive  p  and  a 
positive  T‘  <  T  such  that  for  l$-y0ly<  p,  equation  (2.1)  has  a  unique 
solution  u  6  C°([0,T'];W)  n  C1 ( [0,T' ] ;X).  p  and  T'  depend  only  on  K  and 
the  distance  of  yQ  from  the  boundary  of  W. 

The  second  theorem  stated  in  this  chapter  is  a  continuous  dependence 
result.  It  generalizes  Theorem  II  of  [12]  insofar  as  it  allows  S  to  depend 
on  w.  We  adopt  the  following  assumptions: 

(S2)  There  is  an  open  set  W*  c  z  such  that  W  c  w'  and  the  following 
holds.  The  definition  of  S(t,w)  6  B(Y,X)  can  be  extended  to  w  6  W' . 
Moreover,  we  have  uniformly  on  [0,T]  *  W': 

Is(t.w)  -  <  KIW-»I2  ,  WSIt.if')  -  DwS(t,w)lz>B(r  xJ 

<  Klw’-wl  ,  ID  S(t,w’)  -  D  S( t,w) I  <  Klw’-wl  . 

ti  X.  t  1  I  A  U 


(A5) 
(A6) 
( f  2 ) 


Here  Dw  and  Dt  denote  the  derivatives  w. r.  to  t  and  w. 
I B(t,w’ )  -  B(t,w)lx  <  Klw’-wly 


lA(t,w)  -  A(  t,  w* ) I  <  Klw'-wl 

I  ft*  I 

lf(t,w')  -  f(t,w)ly  <  Klw'-wl^. 


Let  us  now  consider  a  sequence  of  evolution  problems  (n  6  ■) 
(2.2)  un  =  A(t,un)un  +  f(t,un),  0  <  t  <  T,  un(0)  =•  $n. 


Theorem  2.2.: 

Assume  (A1)-(A6),  (fl),  (f2)  are  satisfied  uniformly  in  n  and  assume 
(SI),  (S2).  (The  operator  S  shall  not  depend  on  n. )  Moreover,  assume 
.♦n-y0ly<  P,  **-Yq\  <  P  with  P  as  of  Theorem  2.1.  Finally,  assume  that 
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for  t,w  €  [0,T]  x  W 

An(t,w)  ♦  A( t,  w)  strongly  in  B(Y,Z> 

Bn(t,w)  ♦  B(t,w)  strongly  in  B(X) 

fn(t,w)  ♦  f(t,w)  in  Y 

as  n  ♦  •.  If  $n  ♦  $  in  the  Y-norra  as  n  ♦  •,  then  there  is  a  T"  <  T 
such  that  (2.2)  has  a  solution  un  6  C1 ( [0,T" ] jX)  n  c°( [0,TM ] ; W)  for  any 
n.  Moreover,  un(t)  ♦  u(t)  in  Y,  uniformly  for  t  6  [0,T"),  where  u  is 
a  solution  of  (2.1). 

Proof: 

The  proof  essentially  follows  the  same  line  of  argument  as  that  in  [15]. 
Theorem  2.1  yields  the  existence  of  solutions  to  (2.2)  and  the  limiting  equa¬ 
tion  (2.1)  on  some  interval  [0,T‘ ]  with  T'  independent  of  n.  It  is  more¬ 
over  proved  precisely  as  in  [15]  that  u11  ♦  u  uniformly  in  t  in  the  X- 
norm.  To  prove  convergence  in  the  Y-norm,  we  rely  on  Theorem  IV  of  [14]. 

This  involves  estimating  a  number  of  terms.  Most  of  those  estimates  go  as  in 
[15]  or  are  straightforward,  and  ray  exposition  will  focus  only  on  those  terms 
that  present  difficulties.  As  in  [15],  it  is  used  that  un  solves  the  linear 
equation 

(2.3)  u  *  A  u  +  f  , 

with  An  =•  An(t,un)  and  fn  =  fn(t,un).  The  limit  u  solves  the  linear 
equation 

(2.4)  u  =  Au  +  f 

with  A  =  A(t,u),  f  =  f(t,u).  From  [14],  Theorem  IV,  we  have  the  estimate 
,un-u,a(,Y  <  k( ■  Y  +  <fn-fyY> 

+  K(l(Sn(0)-S(0))$lv  +  I (Sn-S)f«.  v  +  l(Sn-S)ul  v 
+  K(l (Bn-B)sul,  v  +  !(Cn-C)Sul.  v) 

I  j  A  I  f  A 

+  Kd  (Un-U)(5  9  <|>  9  g) I  v. 

9  X 
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V 


Here  Sn  denotes  S(t,un),  Bn  denotes  B(t,un).  The  symbol  C  stand-;  for 
SS  1 .  The  U,Un  are  the  evolution  operators  associated  with  A,An. 

Finally,  denotes  S(0)$,  and  g  stands  for  Sf  +  (C-B)Su.  The  indices 

4  00  ^ 

1  and  00  indicate  the  L  -  and  L  -norms  on  the  interval  tO,T"].  On  the 

right  hand  side  of  (2.5),  the  term  0 4>n-(^ ■  ^  converges  to  zero  by  assumption 

and  we  are  left  with  seven  more  contributions.  Of  these,  the  first,  fifth  and 

seventh  have  been  dealt  with  in  [15],  and  no  change  in  the  argument  is  needed 

here.  The  second,  third  and  fourth  term  are  estimated  in  terms  of 

and  lun-ul  _  by  virtue  of  (SI  ).  Now,  note  that 
Z 

•un-ui  „  <  K(l*“-*l_  +  lun-ul,  )  <  K(l«n-4>l 
®fZ  Z  1  *  Z  x 

+  I Anu°  +  f n  -  Au  -  f  I  ). 

i  ,Z 

The  last  term  will  be  estimated  below.  We  may  thus  focus  on  the  term 
l(Cn-C)Sul,  v  =  KSn(Sn)"1  -  SS_1)Sul,  v 

<  l(Sn-S)ul1x  +  lSn((Sn)*,-S-1 )Sul1  x. 

By  (S2),  S°  is  bounded  in  B(Y,X),  and  by  (SI),  l(Sn)  1 -S  'l  can  be 

X#  x 

estimated  by  lun-ul  This  takes  care  of  the  second  term.  For  the  first 

<",Z 

term,  observe  that 

Sn(t)  -  ~  S(t,un(t))  =  S(t,un(t))  +  D  S(t,un(t))u"(t) 
dt  u 

=  S(t, un(t))  +  D^St t,u°( t) ) (An ( t,un )u°  +  fn(t,un)) 

and  likewise 

S(t)  *  S(t,u)  +  D^S( t, u )  (A(t,u)u  +  f(t,u)). 

We  have 

l(S(t,un)  -  S( t,u) )ul  <  «S(t,un)  -  S( t, u ) I  v#ul 

<  K«un-u«  „  <  KlUn-u«  „ 

1,Z  1,Y 

and 
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r 

I 

t 

I (D  S(t,un(t))  -  D  S(t,u(t) ) ) (An(t,un(t) )un(t)  +  fn(t,un(c) ) )l 

<  .(VU,un(t)>  -  VU'“U,"l.z,B(Y,x)'''V'  *  f"<.  z 

<  U»n-ul,(!,  <  K.an-».1Y. 

Finally, 

'V"'"'”1'"  fn-AU-f  )P5  (B(Y(x)  ‘  'DuSlt'U<t”'-.Z.B(,,X> 

« , n  n  _ n  _» 

•  I A  U  +f  -Au-f 8  . 

1  ,  z 

This  last  term  can  be  estimated  by 

lAn(t,u)u  -  A(t,u)ul  +  Hfn(t,u)  -  f(t,u)l 

1  /  Z  i  f  z 

+  l(An(t,un)  -  An(t,u))ul  +  i An ( t, un ) (un-u) 8 

1  f  Z  1  t  z 

+  lfn(t,un)  -  fn(t,u)l  . 

It  follows  from  Assumptions  (A6),  ( f 2 >  that  the  last  three  terms  can  be  esti¬ 
mated  by  lun-u*1  The  first  two  contributions  converge  to  zero  because  of 


the  assumptions  of  the  theorem.  This  concludes  the  proof 


1 

In  this  chapter,  we  shall  apply  the  preceding  results  to  differential-  j 


delay  equations  of  the  following  form 


(3.1  ) 


U ( n- 1  )  =  ^  f(VUxt'“*'Ux<n-4>>Ux(n-2))> 


+  h(u  ,u  .*»*, u  . 

x  xt  x(n-3)  xx(n-3) 


+  k  ( u  ,  •  •  • ,  u  ,,;u  ,  •  •  • ,  u  . .  i  n )  +  $( t , x ) . 

x  x(n-3)  xx  xx(n-4) 

Here,  the  index  (k)  stands  for  k-fold  differentiation  with  respect  to  the 
time  variable  t,  and  the  hat  denotes  the  past  history:  u*(t,x) (S)  = 

■V  <M  <>✓ 

ux(t+S,X)  for  S  e  (-",0].  The  f,h  and  k  are  smooth  functionals  on  a 
history  space,  the  topology  of  which  will  be  specified  later. 
Differentiating  (3.1)  with  respect  to  x  ,  we  obtain 


(3.2) 


u  ,  .  =  n(f(u  ,u  .•••,u  ,  A.)u  .  _.) 

x(n-1)  x  xt  x(n-4)  x(n-2)  xx 


*  A  “  .  3h 

+  h(u  ,u  -•••,u  ,  ,.;h)u  .  , ,  +  -r—  •  u  , 

x  xt  x(n-3)  xxx (n-3)  ox  xx(n-3) 


.  3d> 

where  ?  =  r1.  In 
dx 


+  (k(u  ,***,u  ,,;u  ,  ••*,u  ,  „,;n))  +  4*  ( t ,  x ) , 

v  x  x(n-3)  xx  xx(n-4)  'x 

the  following,  we  need  only  be  concerned  with 
equations  of  the  form  (3.2).  For  applying  the  results  of  Chapter  2,  it  is  con¬ 
venient  to  rewrite  (3.2)  as  a  system  of  equations.  Let  us  put  v^  =  ) , 

wk  =  uxx(k)’  We  thus  obtain  the  following  system  equivalent  to  (3.2) 
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Jc  -  0 , 1 ,  •  •  •  ,n-3 


vkt  “  vk+1 


wkt  “  wk+1 


0, 1 ,  • • • , n-4 


(3.3) 


A  A  A  A  0, 

/  _  ^  *  n(f (v  -• •  •  f v  Jv  _)  +  h(v_,*##,v  _;n)w  _  +  *r—  •  w 

n-2,t  0  n-4  n-2'xx  0  n-3  n-3,x  3x  n- 

A  A  A  A 

♦  (k(v0'***'Vn-3JV‘“'Wn-4;n)}x+  ♦ 


wn-3,t  =  vn-2, x* 

It  will  be  advantageous  to  make  some  further  substitutions.  Let  us  put 

►  —  •  k.  Then  we  obtain  the  following  £ 
/h 

k  “  0 , 1 , • * • , n-4 


v '  *  fv  w'  B 

f/h  w 

n-2  n-2  n-3 

n- 

V  * 

v 

kt 

k+1 

v  * 

n-3  ,t 

v*  _ 
n-2 

f 

w  * 

w 

kt 

k+1 

w  .  = 

-J-W 

n-4,t 

f/h 

(3.4)  Vn-2,t 

n  f  v'  . 
n-. 

w'  _  * 

(—  w 

f/h  n 
/h  v'  . 

n-3  ,t 

rM 

l 

c 

0 , 1 , • • •  ,n-5 

_  k 
h 

+  /h 


n-3  ,x 


♦  <* '  a  - » g> 


-  5  >  -  *  • » •  <£  >  ♦  ft 

/h  3f  ,  3  , 

•  r  r:  v'  +  r  (f»h)  • 


+  * 


Next,  let  us  define  the  history  spaces,  in  which  (3.4)  will  be  analyzed. 
Definition  3.1; 

For  a  given  Banach  space  Z,  let  H  (Z)  denote  the  space  of  all  func¬ 
tions  (-*,0J  +  Z,  which  are  the  sum  of  a  constant  element  Zg  e  Z  and  a 

function  z(t),  which  is  square  integrable  and  has  a  square  integrable  deriva- 

~2 

tive  (in  the  Bochner  sense).  Analogously,  let  L  (Z)  be  the  space  of  all 
functions  ( — 00 , 0 ]  +  Z,  which  are  the  sum  of  a  constant  and  a  square 
integrable  function. 
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In  particular,  shall  denote  H^(R)  and  shall  denote  H  (H'(R)) 

where  Hm(R)  is  the  Sobolev  space  of  all  functions  R  ♦  R,  which  have  n 

^2  ^2  n 

square  integrable  derivatives.  Analogously,  let  =  L  (H  (R)). 

Remarks ; 

1.  In  [23],  I  used  the  space  (the  space  of  bounded  continuous  func- 

b 

tions  having  a  limit  at  -®.  The  reason  why  this  space  cannot  be  used 
here  is  that  it  is  not  reflexive  as  required  by  Kato's  theory. 

2.  The  choice  of  the  space  seems  to  impose  rather  restrictive  conditions 

on  the  given  history.  However,  as  in  [23],  one  can  allow  histories  in 
more  general  spaces,  e.g.  "fading  memory"  spaces  [7],  by  reducing  the 
problem  to  one  that  has  a  history  in  H1 ,  but  is  equivalent  to  the  given 
one  for  t  >  0.  The  only  modification  necessitated  by  this  is  that 
f,h,k  must  be  allowed  to  depend  explicitly  on  x  and  t.  This 
modification  presents  no  major  difficulties. 

As  in  [23],  we  define  a  shift  operator  T-  on  H,:T„$(t)  *  $(t+s)  for 

s  e  (--,o] . 

Our  assumption  on  f,h,k  in  (3.3)  are  as  follows: 

(i)  The  mappings  fstH^)11  3  ♦  R,  h^H1)11  3  *  R  ♦  R  and 

ktfH1)^11  3  *  R  +  R  are  smooth  (i.e.,  sufficiently  often 
continuously  differentiable)  and  the  induced  operators  f,h,k 

a  a  ^  ^ 

defined  by  f(<p)(S)  =  f(T  $)  map  into  H  and  depend  again  smoothly 

on  their  arguments.  Moreover,  f  and  h  take  strictly  positive 

values:  f>e>0,h>e>0.  k  vanishes  if  its  arguments  are 

zero.  Moreover,  the  Frechet-derivative  Df  is  a  linear  operator 

from  (L  )  into  R,  which  depends  smoothly  on  the  arguments  of 

f  (in  the  topology  of  H  ),  and  the  corresponding  operator  Df  maps 
~2 

into  L  and  again  depends  smoothly  on  the  arguments  of  f. 

Analogous  conditions  hold  for  h  and  k. 


The  following  lemma  is  easily  proved. 

Lemma  3.2: 

If  (i)  holds,  with  a  sufficient  degree  of  differentiability,  then  f 

defines  in  a  natural  way  (acting  pointwise  in  the  space  variable  x)  a  smooth 

operator  from  ( H 1 ) n  3  into  H1  +  R.  Here  m  is  a  given  integer  greater  or 
mm 

equal  to  1.  The  same  holds  for  h,k,  and  an  analogous  statement  also  holds 

*  ~2  n-3  ~2  2n-3  ~2 

for  Df,  Dh,  Dk  (regarded  as  linear  operators  (L  )  or  (L  )  ♦  L  ) . 

Remark: 


Since  we  are  concerned  with  existence  theorems  locally  in  time,  it  is 
clearly  sufficient  that  condition  (i)  holds  in  a  neighborhood  of  the  pre¬ 
scribed  initial  condition. 


3h  p 

Let  us  also  note  that  -r—  =  l  D.h(v  , 


n-2 


i=1 


',v  ,;h)w,  , ,  where  D. 

n-3  1-1  l 


denotes  the  Frechet  derivative  w.r.  to  the  ith  argument.  Analogous 
manipulations  are  possible  for  the  other  x-  and  t  derivatives  of  f,h 
and  k  that  occur  in  (3.4).  As  in  [23],  let  us  assume  that  the  given  initial 


history  up  to 

time  t=0 

satisfies 

the  equations  (this  can  be  achieved 

appropriately 

changing 

♦.  We  can 

then  rewrite  equation  (3.4)  in  the 

following  abstract  form 

V. 

=  V, 

k  =  0,1, • 

• • ,n-4 

kt 

k+1 

w 

kt 

=  Wk+1 

k  =  0,1,* 

•*,n-5 

(3.5) 


n-2 


n-3,t 


w 


n-4,t 


n-2  ,t 


n-3 

sz 

At  A 

n  f  v* 


k 

h 


n-2, xx 


+  /h  w 


n-3,  x 


+  •  •  •  + 


n-3  ,t 


■  Ay- 


n-2,  x 


+  •  • « 
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where  f,h  etc.  are  defined  as  in  Lemma  3.2.  Now  let  the  spaces  X,Y  be  as 

follows:  X  =  (h]  n  L2)2n  Y  =  (H^  n  l2  Yn  ^  ,  where  m  is  an  odd  number 
14  m  m+J 

~1  ~2  2n-5 

greater  or  equal  to  3.  Finally,  let  Z  =  (H  „  ft  L  .)  .  When  identifying 

ra-2  1 


(3.5)  with  the  abstract  system  (2.1),  we  incorporate  in  A  only  those  terms 

that  contain  derivatives  w.r.  to  x,  everything  else  is  included  in  f.  With 

this  identification,  the  conditions  (A3),  (A4),  (fl),  (A6)  and  (f2)  are  rather 

<>*2 

obvious,  (note  that,  for  m  >  3,  H  n  L  _  is  a  Banach  algebra)  provided 

m  m+  j 

that  the  following  holds: 

(ii)  ♦  takes  values  in  H  (it  follows  automatically  that  it  is  continuous 

on 

into  H  ),  and  the  initial  condition  at  t=0  lies  in  Y. 
m 

For  verifying  the  remaining  conditions,  we  have  to  study  the  operator 


A(f,h),  defined  by 


A(f,h)(v,w)  =  (n  f 


v  +  Jh  w  ,  v  ). 


/  o  »  2 

The  operator  S  is  defined  by  S  ■  (,f  — -  -  A j  for  A  e  R  large 

3x 

enough.  With  this  choice  of  S,  conditions  (Si),  (S2)  are  obvious. 
Moreover,  one  sees  easily  that 

.  .2  A  A  AAA  ^2  A  A 

((f  — T  -  A]A(f,h)  -  A( f ,h ) (f  —  -  A) )(v,w) 

3x  3x 

yields  an  expression  involving  only  first  and  second  order  derivatives  of 

AAA  A 

v,w,f  and  h.  From  this  it  is  not  difficult  to  conclude  (A2)  and  (A5), 


(Note  that  TnA-ATn  =  Tn_1 (TA-AT)  +  Tn-2(TA-AT)T  + 


and  apply  this 


*  J*  A  A  A  A 

with  T  =  (f  — r  -  Aj).  For  (Al),  we  have  to  show  that  Re (A(v,w ) , ( v,w) )  „ 

3x2  (kY 

A  A  A  A  1 

C( (v,w) , ( v,w) )  which  follows  from  a  simple  integration  by  parts  in  the 


x- variable. 
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We  have  thus  proved: 

Theorem  3.2. 

Let  (i),  (ii)  be  satisfied.  Then  there  is  a  T  >  0  such  that  (3.5)  has 

"  '  *  1  0  * 
a  solution  U  -  (vA,v  ,•••,«'  -)  e  C  ([0,T)>X)  n  c  (I0,T];Y).  U  depends 
u  i  n-^ 

continuously  on  he  [0,h0]  in  the  norm  of  Y. 
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Stretching  of  Filaments  of  Viscoelastic  Liquids 


We  are  studying  the  motion  of  an  infinitely  extended  filament  of  an  incom¬ 
pressible  viscoelastic  liquid  under  the  influence  of  a  longitudinal  body 
force.  It  was  shown  in  [233  that,  if  the  filament  is  thin,  this  problem  can 
be  modelled  by  a  one-dimensional  approximation,  where  only  longitudinal  mo¬ 
tions  need  to  be  studied.  Let  u<x,t)  denote  the  position  of  a  fluid  par¬ 
ticle  at  time  t,  which  is  at  the  position  x  in  certain  reference  state. 

For  simplicity,  this  reference  configuration  is  chosen  to  be  one  in  which  the 
filament  has  uniform  thickness.  I  showed  in  [23]  that  the  evolution  of  u  is 
governed  by  the  following  equation 

(4.1)  pu^  =  I—  (u  x1 1  -  u  2x22)  +  $. 
tt  ox  x  x 

In  this  equation  P  denotes  the  density  of  the  fluid  (i.e.,  a  constant),  4 

11  22 

is  the  given  body  force,  and  it  ,it  are  the  longitudinal  and  transverse 
components  of  the  convected  extra  stress  tensor  (i.e.,  not  including  the  pres¬ 
sure,  which  was  eliminated  in  the  derivation  of  (4.1)).  The  tensor  it  is 
related  by  a  constitutive  law  to  the  right  Cauchy-Green  tensor  y  (in  our  no¬ 
tation  we  follow  [22]).  In  the  approximation  leading  to  (4.1),  y  is  given 
by 


In  the  following,  we  discuss  various  constitutive  laws  that  have  been  sug¬ 
gested  in  the  rheology  literature  and  the  corresponding  equations  (4.1)  that 
they  lead  to.  tt  will  be  shown  that  all  these  equations  can  be  transformed  to 
the  form  (3.1).  In  particular,  we  shall  check  the  positivity  of  the  func¬ 
tions  f  and  h.  (It  will  always  be  understood  that  u  is  the  sum  of  a 
given  function  uQ(x,t)  and  a  function  tending  to  zero  appropriately  as 


x  ♦  if,  moreover,  ux  is  always  assumed  uniformly  positive.  In  comparison 
to  Chapter  3,  the  variable  called  u  there  will  be  identified  with  u- 
uQ(x,t)  in  this  chapter.)  The  notation  used  in  the  original  papers  cited 
here  is  often  different  from  ours,  and  we  have  transcribed  the  constitutive 
laws  appropriately.  Tables  of  some  (but  not  all)  constitutive  assumptions 
discussed  here  can  be  found  in  [2]  and  [22). 

a)  The  rubberlike  liquid  of  Green  and  Tobol sky  [10)  and  Lodge  [T9],  [20]  and 
modifications  of  Ward  and  Jenkins  [27)  and  Lodge  [21]. 

In  these  theories,  the  constitutive  law  has  the  following  form 

*  *  ~n  Tt  *  +  t_s (s)ds  -  /t  b ( t-s ) Y  1 (t)Y(s)Y  1 (t)ds 

^  _oo  _00 

*  ft  c(t-s)Y-1  <s)Y(t)Y  1  (s)ds  +  ft  d(t-s)  (Y(t):  Y  NsJJY  ^sjds. 

_00  —49 

The  first  term  is  a  Newtonian  contribution,  the  second  is  the  one  given  by  the 
rubberlike  liquid  theory  [10],  [19],  [20).  The  third  term  accounts  for  a 
modification  suggested  by  Ward  and  Jenkins  [27],  and  the  last  two  represent 
corrections  of  Lodge  [21]  (Lodge  finds  c=2d  from  a  molecular  theory,  but  we 
shall  make  no  use  of  this).  With  this  constitutive  law,  Equation  (4.1)  as¬ 
sumes  the  following  form 

2 

Pu  *  3h  a~5T  ( - — )  +  I— {u3(t)  •  /t(c(t-s)+d(t-s))u  4(s)ds 

tt  3x9 t  ^  u  '  9x  x  '  x 

x  *" 

+  u  (t)  •  /fc  a(t-s)u  2(s)ds  +  /t  b(t-s)+d(t-s))u  1(s)ds 

X  _oo  X  _oo  X 

-  u  2(t)  •  Jt  a(t-s)u  (s)ds  -  u  3(t)  •  /*"  (b( t-s )+c( t-s )+2d( t-s ) )u2(s )ds} 

X  X  X  ^  X 

jto  »oo 

+  ♦  . 

In  order  to  obtain  the  form  (3.1),  we  differentiate  this  once  with  respect  to 
time.  This  yields 
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1  -  xx  ,  xt  ^2 

PU  =  3n  — —  •  u  .  .  -  6n  — r—  u  +  u  .  •  1-1 20——  +  3u  • 
ttt  2  xxtt  3  xtt  xxt  3  x 

u  u  u 

XX  X 

/fc  (c(t-s)  +  d<t-s))u  4  (s)ds  +  /*  a(t-s)u  2(s)ds  +  2u  3  • 

_•  x  _•»  x  x 

/ta(t-s)u  (s)ds+3u  4(t)  •  (b( t-s ) +c ( t-s )+2d( t-s ) )u2(s)ds}  + 

»  V  tf  *  V 


Here,  as  always  in  the  following,  the  dots  indicate  terras  involving  only  lower 
order  derivatives  of  u  and  the  derivatives  of  ♦  (♦  is  always  assumed 
“smooth  enough").  We  have  assumed  here  that  the  kernels  have  derivatives  in 
L1  so  that 

ft  a(t-s)f(s)ds  *  /*  a  '.t-s  )f  (s  )ds  +  a(0)f(t). 

The  equation  above  clearly  has  the  form  {3.1),  and  the  coefficient  of  uxxt 

is  positive  if  the  kernels  are  positive  and  0  is  small  enough. 

b)  The  model  of  Kaye  [17)  and  Jsernstein,  Kearsley  and  Zapas  [1]. 

In  this  model,  the  constitutive  law  has  the  form 

*  -  Jt  a(t-s,It  »I2)T~1  (s)ds  -  J*  b(t-s,I},I2)Y”1(t)Y(s)ir‘1  (t)ds. 

.00  .00 

I1  and  I2  are  the  invariants  of  Y  1(t)Y(s):I1  *  tr(Y(t)Y  Ns))  and  i2  « 

-1  2  -2 
tr(Y  (t)Y(s)).  In  our  special  problem,  we  thus  have  I1  =  u^ttju^  (s)  + 

-1  2  -2  -1 
2u  (t)u  (s)  and  I„  »  u  (s)u  (t)  +  2u  (s)u  (t).  Both  are  thus  functions 
xx  2  x  x  xx 

ux(t) 

of  the  single  variable  I  *  — ; — r  «  and  we  shall  use  the  obvious  notation 

ux(s) 

a (t-s, I),  b(t-s,I).  The  dynamic  equation  (4.1)  assumes  the  form 

Pu..  *  ■»—  (u  (t)  ft  a(t-s,I)u  2(s)ds  +  ft  b(t-s,I)u  Nsjds 

tt  OX  X  —  X  X 

.00  .00 

-u  2(t)  ft  a(t-s,I)u  (s)ds  -  u  3(t)  •  /t  b(t-s,I)u2(s)ds)  + 


By  differentiation  with  respect  to  time,  we  find 
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Pu.  »  u  .{/*  a(t-s,I)u~2(s)ds  4-  2u~3(t)  ft  a(t-s,I)u  (s)ds 

ttt  xxt  '  X  X  '  X 

— »«>  _oo 

+  3u  4  /*  b(t-s,I)u2(s)ds  +•  u  (t)  /t  (t-s,I)u  3(s)ds 

X  X  x  w  OX  X 

+  /*  (t-s,I)u~2(s)ds  -  U;2(t)  / t  Iy  (t-s,I)ds 

-  u“3(t)  /*  (t-s,I)ux(s)ds}  +  •••  . 


A  Newtonian  term  can  be  added  to  this  as  before.  Suppose  the  kernels  a,b 

are  positive.  Then  the  coefficient  of  uxxt  is  positive  in  two  cases: 

3a  3b 

a).  If  ,  yj  are  small,  i.e.,  if  the  model  is  considered  a 

perturbation  of  the  Hard -Jenkins  model. 

3a  3a  3b  3b 

0)  If  Ti~'  TT~  are  Posi-t*ve*  It  be  interesting  if 

this  condition  has  a  physical  interpretation, 
c)  The  Bird-Carreau  model  [31,  [5], 

In  this  model,  we  have 

*  -  <1  +^)  /ta(t-s(I(s))Y'1(s)ds  -  J*  a(t-s,I(s))Y"1(t)Y(s)Y"1(t)ds 
^  -*•  ^  _• 

2  /  V 

•  -1  •  -1  Uxt(8> 

where  I(s)  *  tr(Y(s)Y  (s)Y(s)Y  (s))  *  6 


Ux(8) 


This  leads  to  the  equation 


Pu  *  (1  +  |—(u  Jt  a(t-s,I(s))u  2(s )ds 

tt  2  ox  x  '  x 


-  u  2 ( t )  /fc  a(t-s,I(s))u  (s )ds} 

X  J,  X 

+  \  Jt  a(t-s,I(s) (u^1 (s)ds 

-  u  3(t)  Jt  a(t-s,I(s) )u2(s)ds}  + 

X  X 
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I 


Differentiating  this  w.r.  to  time,  we  find 


Puttt  =  uxxt t < 1  + -|){  /ta(t-s,I(s)u^2(s)ds  +  2uw  (t)  Jt  a( t-s , I ( s ) )u^(s)ds} 


+  j  •  3ux4  Jt  a (t-s, I ( s  )  )u2 (s )ds) 


+  lO  +  ^){u  ( t)  jZ  4r-  ( t-s,  I(s  ))I  (s)u  2(s)ds 


3t 


2  x 


31 


u2  ft  It-  ( t-s,I(s)  )I  (s)u  (s)ds}  +  4  {  /fc  4f-(  t-s,  I(s ) ) 


X  X 


e  ,  ft  3a, 


I  (sJu'Vjds  -  u-3  fZ  |4  ( t— s ,  I  ( s ) )  I  (s  )u2(s  )ds}]  + 

XX  X  ^  <»I  XX 


A  second  differentiation  yields  (the  kernels  are  assumed  to  be  twice 
differentiable  w.r.  to  time 

Putttt  =  uxxttl(1  +f){  I*  a(t-s,I(s))u“2(s)ds 
+  2u  3  ft  a(t-s,I(s))u  <s)ds} 

x  __  X 

+  4  3u  4  Jt  a(  t-s  ,I(s)  )u2(s)ds]  +  •••  . 

x  X  ^  X 

For  a  positive  kernel,  the  coefficient  of  uxxtfc  is  positive. 

d)  The  Carreau  Model  B  [41. 

In  this  model,  it  is  assumed  that 

“  (1  -j)  ft  exp(-  ft  f(I(r))dr)Y  ^sjds 

s 

“  j  !t  exp(-  ft  f (I(r) )ds)Y_1 (t)Y(s)Y_1 (t)ds, 
s 

where  I  has  the  same  meaning  as  in  the  Bird-Carreau  model.  We  thus  obtain 
the  following  equation 

Putt  *  JZ  {  ft  exP<-  /*  f(I(r))dr)[(1  +  f)(u  (t)u"2(s)  -  u~2(t)u  .(s)) 

“  X  X  ° 


X  X 


+  4<u  ’(s)  -  u  3( t)u2(s ) )ds ) }  +  $. 
2  X  XX 
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The  integral  converges,  if  f  takes  strictly  positive  values.  Differenti¬ 
ating  with  respect  to  time,  we  obtain 

Pu.k.  =  u  .{  ftex  p(-  /fc  f  {I(r )  )dr)  [  (1  +  zr)(u  2  ( s )  +  2u  3(t)u  (s)) 
ttt  xxt  L  2  x  xx 

+  4  3u-4(t)u2(s) ]ds}  +  4— {-  /t  exp(-  / t  f(I(r))dr) 

2  x  x  at  '  ' 

s 

•  /t  f'(I(r))I  (r)dr [( 1  +  ~)(\i  (t)u-2(s)  -  u“2(t)u  (s))  + 
x  2  x  x  xx 

s 

+  -^(u  1 (s)  -  u  3(t)u2(s ) ) Ids}  +  ••• 

2  x  xx 

The  second  differentiating  with  respect  to  time  yields 

Pu  *  u  {  /fcex p(-  Jt  f (I(r) )dr  1(1  +~)(u  2(s)  +  2u  3(t)u  (s)) 
tttt  xxtt  _  2  x  xx 


+  ~  3ux4(t)u2(s ) Ids  -  12 


u  (t) 

XC  -  /  exp(-/  f(I(r))dr)f*{I(t)) 


u  (t)  -®  s 

X 

•  [(1  +  4)(u  (t)u~2(s)  -  u“2(t)u  (s)) 

*  X  X  x  x 

+  ^-(ux,(s)  -  u^3( t )u2(s ) )  Ids}  +  ••• 

The  coefficient  of  uxxtt  *-s  positive  under  the  restriction  that  f'(l)/T  is 
not  too  big. 

e).  The  Leonov  model  [18], 

This  model  does  not  explicitly  give  the  stress  as  a  functional  of  the 
strain  history.  Instead  it  is  given  by  a  system  of  equations  as  follows 


*  ’  l  "!’‘’‘,1k'V  %  - 


7t*ck  *  “  ‘fxtItk'I2k,tcX  Yck  "  3  Xlkck  1 


’9X(Itk,I2k)<3  X2k  ck  "  y  1 
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where  I1k  =*  tr(c”  Y),  I2Jc  =  tr(Y~  c^).  The  tensors  ck  satisfy  the  restric¬ 
tion  det  =  1  (it  can  be  shown  that  det  is  an  inv  .iant  of  the  evolu¬ 

tion  equation).  In  Leonov's  paper,  the  analogue  of  c^1  xs  called  ck,  we 

(It )  (k ) 

have  changed  this  for  consistency  of  notation.  The  ,W  ,  w'*tc'^k  are 

positive  scalar  functions,  they  are  not  independent  in  Leonov's  model.  It  is 
convenient  to  introduce  d^  =  c^  Y*  With  this,  the  constitutive  equation 


becomes 


•  ■  l  -  W21'>‘I1R'12X,dir,T 

k 

-n  "“n-hi* 


-1.-1 


(4.2) 


“  ■fktIlk'12kJtdk  "  3  I1kdkt 


•Vhk'VTVk  -  "  *  V'’; 


I1k  and  I2k  are  the  first  and  second  invariants  of  dk,  and  we  have  det  d^ 

»  1 .  If  fk,  gk  have  positive  values,  then,  for  y  =  0,  the  solution  d^  = 

id  is  an  exponentially  asymptotically  stable  solution  of  (4.2)  when  this 

—  1  * 

equation  is  restricted  to  {d^ldetd^  =  l}.  Consequently,  if  Y  Y  +  0  as 
t  -*■  then  on  some  interval  (-®,T)  there  is  a  unique  solution  dk  which 

converges  to  the  identity  as  t  +  Whether  this  solution  can  be  con¬ 

tinued  up  to  t  =  0  depends  on  the  form  of  f^,  gk  and  the  history  of  Y. 

We  shall  assume  that  (4.2)  has  a  solution  up  to  t  =  0.  Then  this  solution  is 
a  smooth  functional  of  the  histories  of  Y  and  Y:  dk  =  (Y,Y)»  From 

(4.2)  or,  resp. ,  its  differentiated  version,  we  also  find  functional  relation- 

A  A 

*  *  .  *  •  -1,, 

ships  of  the  form  d^  =  F2(Y,Y),  -  F3(Y,Y)  +  c^Y  Y.  For  the  filament 

problem,  Y  is  a  diagonal  matrix,  and  so  is  d^.  Let  us  denote  the  11-  and 
22-components  of  d^  by  dk1 ,  dk2.  The  dynamic  equation  (4.1)  reads  now  as 
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follows 


Putt  W!k)(111c'I2lc)dk1  “  W2*/(I1k'I2k,\i 

k 

-  W,  lI1k'I2k)dk2  +  W2  (ITk'I2k,dk2I 

+  3n  w(In'I2i)  IT  l--T)}  + 

x 

Retaining  only  terms  of  the  highest  differentiation  orders,  we  obtain  by  two¬ 
fold  differentiation 

"“tttt  ’  3"  ""n'1!.1  -T  “xxttt  *  Vtt  0(0>  *  'Wt(0<n) 

u 

*  *  3  “2k>d*J  *  a,!k>d«»  *  -  dk2|2 

*  -  C’2  *  2l»,lk’  *  »"»dk,  -  \2>ld«  • 


.(k), 


-1 


3w 


(k) 


(K )  1 

Here  W.  .  stands  for  yj -  .  For  small  h,  the  coefficient  of  uxxfct  is 

13  jk 

•  ..  ,  ._  u(k)  „(k)  u(k)  (k)  ...  ,  ,„(k)  .  „(k).2 

positive  in  particular  if  ,W^  ,W^  'W22  are  Posltlve  and  (w12  +  W21 

(k)  (k) 

<  4  W22  .  This  corresponds  to  inequality  [1.33]  in  Leonov's  paper, 

f ).  The  models  of  Johnson  and  Segalman  [13]  and  Chang,  Bloch  and  Tschoegl  [28]. 
This  model  is  described  by  the  following  system 

it  =  ft  a(t-s)G(s,t)Y  1  <s  )GT(s,  t)ds 

-00 

yy  *  -cpf  ( t )  Y  ( t )  G 


G(t,t)  =  id 

«  ot  G  Y-1(s)Y<s). 

9s 

The  parameter  a  ranges  between  0  and  V2  •  For  our  problem,  Y  is  diag¬ 
onal,  whence  Y{t'),Y(t")  commute  for  any  t',t".  The  equations  for  G  can 
therefore  be  solved  as  follows. 
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G(s,t)  =  exp(-a  /C  Y  (r)Y(r)dr)  =  Y  (s)Y  (t). 
s 

This  leads  to  an  equation  very  similar  to  the  ones  studied  in  part  a,  and  the 
discussion  follows  closely  the  one  given  there.  We  leave  the  details  to  the 
reader.  For  a  in  the  range  (0,  V2  the  coefficient  of  uxxtt  turns  out 
to  be  positive.  If  a  is  allowed  bigger  than  V2  *  the  type  of  the  equation 
may  change  from  hyperbolic  to  elliptic. 

The  model  of  Chang,  Bloch  and  Tschoegl  is,  for  this  particular  problem, 
equivalent  to  that  of  Johnson  and  Segalman. 


g)  The  model  of  Curtiss  and  Bird  [8], 


This  model  proposes  the  following  constitutive  law 

T 

it  =  -ft  f  a(t-s)[1  +  vT(Y(s  )-Y(t)  )vj  3//2  -  - 

—Oo  ft  /  T  2 

Y (t)  /v  Y  (t)v 

-  H  /t  /  b( t-s ) [ 1  +  vT (Y (s )-Y ( t) ) v]  3^2vTY(t)v  *  — 


/  T  2,  . 
'v  Y  ( t ) v 


Here  is  the  set  Q  ^  55  {v|vTY(t)v  =  l}.  With  £1  denoting  the  unit 


sphere,  this  yields  for  our  problem 


/t  /  a  ( t-s )  [  1  +  (u  2(t)u2(s)-1  )w2  +  (u  (t)u  ^sj-l) 

*  4  XX  1  X  X 

-00  £) 


,  2  2. ,-3/2  2  -2.  .  . 

(w  +  w-) J  '  w,u  (t)dw 
2  3  lx 


-  n  /  b(t-s)[l  +  (u  2(t)u2(s)-1)w2  +  (u  (t)u  1 { s ) — 1 ) 

-  U  XX  1  X  X 

_oo  51 


,  2  2, ,-3/2  2  -3,  •  .  2  2 

W2+W3)  W1Ux  <t>ux<t)(2wrw2-w3)dw. 


22  2-22 
For  *  ,  we  have  the  same  exprssion  with  w^^  replaced  by  W2UX*  when 

inserting  this  into  the  dynamic  equation  (4.1),  we  can  again  achieve  the  form 

(3.1)  by  differentiating  with  respect  to  time.  The  term  involving  u  tt  has 

a  positive  coefficient  proportional  to  n,  the  coefficient  of  uxxt  is,  up 

to  terms  of  0(n) 
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/t  /  a(t-s)[***]  3^2 (w2-w2 )dw  u  2(t) 

-•a  1  x 

♦  4  /fc  /  a(t-s)[***l  5//2(-2u  4(t)u2(8)w2  +  u  \t)u  1(s)(w2+«2)  • 

2  jj  XXiXX^J 

2  2 

•  (w2~w1 )dw. 

It  can  easily  be  checked  that  this  coefficient  is  positive  in  a  neighborhood 
of  the  rest  state. 

Remark: 

When  we  differentiated  equations  with  respect  to  time,  we  have  always 

assumed  that  the  integral  kernels  were  sufficiently  smooth.  Some  of  the 

kernels  suggested  in  the  literature  have  singularities  at  t  =  0  (see  e.g. 

2 

V  “01  t 

( 8] ,  where  a(t)  =  2.  e  )•  A  mathematical  theory  accomodating  such 

a  odd 

kernels  would  be  of  interest.  Experimental  data  on  polymer  melts  (see  e.g. 
[28])  also  seem  to  suggest  that  the  integral  kernel  may  be  singular  at  t  =  0. 


-25- 


REFERENCES 


1.  B.  Bernstein,  E.  A.  Kearsley  and  L.  J.  Zapas,  A  Study  of  Stress 
Relaxation  with  Finite  Strain,  Trans.  Soc.  Rheol.  7  (1963),  391-410. 

2.  R.  B.  Bird,  et.al..  Dynamics  of  Polymeric  Liquids  (2  vol.),  J. 
Wiley,  New  York  1977. 

3.  R.  B.  Bird  and  P.  J.  Carreau,  Chem.  Eng.  Sci.  23  (1968),  427-434. 

4.  P.  J.  Carreau,  Rheological  Equations  from  Molecular  Network 
Theories,  Trans.  Soc.  Rheol.  16  (1972),  99-128. 

5.  I.  J.  Chen  and  D.  C.  Bogue,  Time-Dependent  Stress  in  Polymer  Melts 
and  Review  of  Viscoelastic  Theory,  Trans.  Soc.  Rheol.  16  (1972), 
59-78. 

6.  B.  D.  Coleman,  M.  E.  Gurtin  and  J.  Herrera,  Waves  in  Materials  with 
Memory,  Arch.  Rat.  Mech.  Anal.  19  (1965),  1-19  and  239-298. 

7.  B.  D.  Coleman  and  W.  Noll,  An  Approximation  Theorem  for  Functionals, 
with  Applications  in  Continuum  Mechanics,  Arch.  Rat.  Mech.  Anal.  6 
(1960),  355-370. 

8.  C.  F.  Curtiss  and  R.  B.  Bird,  A  Kinetic  Theory  for  Polymer  Melts,  J. 
Chem.  Phys.  74  (1981),  2016-2033. 

9.  C.  M.  Dafermos  and  J.  A.  Nohel,  A  Nonlinear  Hyperbolic  Volterra 
Equation  in  Viscoelasticity,  MRC  TSRI2095,  Univ.  of  Wisconsin-Madison 
1980. 

10.  M.  S.  Green  and  A.  V.  Tobolsky,  A  New  Approach  to  the  Theory  of 
Relaxing  Polymeric  Media,  J.  Chem.  Phys.  14  (1946),  80-100. 

11.  J.  K.  Hale,  Theory  of  Functional  Differential  Equations,  Springer, 
Berlin-Heidelberg-New  York  1977. 


-26- 


12.  T.  J.  R.  Hughes,  T.  Kato  and  J.  E.  Marsden,  Well-posed  Quasi-linear 
Second-order  Hyperbolic  Systems  with  Applications  to  Nonlinear 
Elastodynamics  and  General  Relativity,  Arch.  Rat.  Mech.  Anal.  63 
(1976),  273-284. 

13.  M.  W.  Johnson  and  D.  Segalman,  A  Model  for  Viscoelastic  Fluid 
Behaviour  which  Allows  Non-affine  Deformation,  J.  Non-Newtonian  Fluid 
Mech.  2  (1977),  255-270. 

14.  T.  Kato,  Linear  Evolution  Equations  of  "Hyperbolic"  Type  II,  J. 

Math.  Soc.  Japan  25  (1973),  648-666. 

15.  T.  Kato,  Quasi-linear  Equations  of  Evolution  with  Application  to 
Partial  Differential  Equations,  in:  Spectral  Theory  of  Differential 
Equations,  w.  N.  Everitt  (ed.).  Springer  Lectures  Notes  in  Math.  448 
(1975),  25-70. 

16.  T.  Kato,  Linear  and  Quasilinear  Equations  of  Evolution  of  Hyperbolic 
Type,  C.  I.  M.  E.  1976. 

17.  A.  Kaye,  Co  A  Note  No.  134,  The  College  of  Aeronautics,  Cranfield, 
Bletchley,  England,  1962. 

18.  A.  I.  Leonov,  Nonequilibrium  Thermodynamics  and  Rheology  of 
Viscoelastic  Polymer  Media,  Rheol.  Acta.  15  (1976),  85-98. 

19.  A.  S.  Lodge,  Proc.  Int.  Congr.  Rheol.  2nd,  1953  (V.  G.  W.  Harrison, 
ed.)  229,  Butterworth,  London. 

20.  A.  S.  Lodge,  A  Network  Theory  of  Flow  Birefringence  and  Stress  in 
Concentrated  Polymer  Solutions,  Trans.  Faraday  Soc.  52  (1956),  120- 
130. 

21.  A.  S.  Lodge,  Constitutive  Equations  from  Molecular  Network  Theories 
for  Polymer  Solutions,  Rheol.  Acta  7  (1968),  379-392. 


-27- 


22. 

23. 

24. 

25. 

26. 

27. 

28. 

29. 

30. 

31. 


MR/db 


A.  S.  Lodge,  Body  Tensor  Fields  in  Continuum  Mechanics,  i.c 

Press,  New  York  -  San  Francisco  -  London  1974. 

M.  Renardy,  A  Semigroup  Approach  to  Equations  with  Infinite  Delay 
and  Applications  to  a  Problem  in  Viscoelasticity,  MRC  TSR#2287,  univ. 
of  Wisconsin-Madison  1981  (to  appear  in  J.  Diff.  Eq. ) . 

J.  C.  Saut  and  D.  D.  Joseph,  Fading  Memory,  in  preparation. 

P.  E.  Sobolevskii,  Equations  of  Parabolic  Type  in  a  Banach  Space, 
Amer.  Math.  Soc.  Transl.  49  (1966),  1-62. 

C.  A.  Truesdell  and  W.  Noll,  The  Non-linear  Field  Theories  of 
Mechanics,  Handbuch  der  Physik  III/3,  Springer,  Berlin-Heidelberg-New 
York  1965. 

A.  F.  K.  Hard  and  G.  M.  Jenkins,  Normal  Thrust  in  Dynamic  Torsion 
for  Rubberlike  Materials,  Rheol.  Acta  1  (1958),  110-114. 

W.  V.  Chang,  R.  Bloch  and  V.  W.  Tschoegl,  On  the  Theory  of 
Viscoelastic  Behaviour  of  Soft  Polymers  in  Moderately  Large 
Deformations,  Rheol.  Acta  15  (1976),  367-378. 

H.  M.  Laun,  Description  of  the  Non-linear  Shear  Behaviour  of  a  Low 
Density  Polyethylene  Melt  by  Means  of  an  Experimentally  Determined 
Strain  Dependent  Memory  Function,  Rheol.  Acta  17  (1978),  1-15. 

M.  Slemrod,  A  Hereditary  Partial  Differential  Equation  with 
Applications  in  the  Theory  of  Simple  Fluids,  Arch.  Rat.  Mech.  Anal. 
62  (1976),  303-321. 

E.  F.  Infante  and  J.  A.  Walker,  A  Stability  Investigation  for  an 
Incompressible  Simple  Fluid  with  Fading  Memory,  Arch.  Rat.  Mech. 

Anal.  72  (1980),  203-218. 


-28 


SECURITY  CLASSIFICATION  OF  THIS  RAGE  (Whtn  Datt  Enfere d) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.  RECIPIENT'S  CATALOG  NUMBER 

O  M  A 

'  J  Jr  J  i 


«.  TITLE  (mxd  Submit)  s.  TYPE  OF  REPORT  b  PERIOD  COVERED 

singularly  perturbed  hyperbolic  evolution  problems  Summary  Report  -no  specific 

WITH  INFINITE  DELAY  AND  AN  APPLICATION  TO  POLYMER  reporting  period _ 

RHEOLOGY  *■  PERF°RMIMG  oRG-  report  number 


7.  author (t) 

Michael  Renardy 


8.  CONTRACT  OR  GRANT  NUMBERf.) 


DAAG29-  80-C-00  41 


».  PERFORMING  ORGANIZATION  NAME  AND  AOORESS 

Mathematics  Research  Center,  University  of 
610  Walnut  Street  Wisconsin 

Madison.  Wisconsin  53706 


It.  CONTROLLING  OFFICE  NAME  AND  AOORESS 

U.  S.  Army  Research  Office 
P.O.  Box  12211 

Research  Triangle  Park.  North  Carolina  27709 


I 


MONITORING  AGENCY  NAME  A  ADDRESSfJf  different  from  Controlling  Office) 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  A  WORK  UNIT  NUMBERS 


1  -  Applied  Analysis 


12.  REPORT  DATE 

May  1982 


IS.  NUMBER  OF  PAGES 

28 


IS.  SECURITY  CLASS.  fof  thlt  report) 


UNCLASSIFIED 


ISA.  DECLASSIFICATION/DOWNGRADING 
SCHEDULE 


l«.  DISTRIBUTION  STATEMENT  (at  thlt  Report) 

Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  the  mbmtrmct  mntered  In  Block  30,  It  different  from  Rmport) 


19.  KEY  WORDS  (Continue  on  reveree  »idm  if  neceeemry  end  Identity  by  block  number) 


Quasilinear  Hyperbolic  Equations,  Differential  Delay  Equations,  Semigroups, 
Singular  Perturbations,  Viscoelastic  Liquids. 


20.  ABSTRACT  (Continue  on  rereree  ride  If  neceeemry  and  Identify  by  block  number) 

•Av'We  prove  an  existence  theorem  locally  in  time  for  quasilinear  hyperbolic 
equations,  in  which  the  coefficients  are  allowed  to  depend  on  the  history  of 
the  dependent  variable.  Singular  perturbations,  which  change  the  type  of  the 
equation  to  parabolic,  are  included,  and  continuous  dependence  of  the  solu¬ 
tions  on  the  perturbation  parameter  is  shown.  It  is  demonstrated  that,  for  a 
substantial  number  of  constitutive  models  suggested  in  the  literature,  the 
stretching  of  filaments  of  polymeric  liquids  is  described  by  equations  of  the 
kind  under  study  here. 


00  i  j  an*7j  1473 


COITION  OF 


65  IS  OBSOLETE 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Dele  Entered) 


